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Abstract. We study the significance of ^Li abundances measured in metal-poor 
halo stars. We explore possible depletion factors for ^Li, defined as the ratio of 
the proto-stellar to the observed abundance, in the three stars where it has been 
detected; to this end, we assume that ^Li/H scales as ^^0/H throughout the galac- 
tic evolution. This assumption is motivated by the recent observations of a similar 
scaling for ^Be, and by the fact that, apart from a — a fusion creation of ^Li, both 
elements are only created in p,a-C,N,0 spallation reactions. We examine possi- 
ble uncertainties attached to the observations and to the modeling of ^Li galactic 
evolution; notably, we include a recent evaluation of the primordial production of 
^Li. The depletion factor Dq in the hottest turn-off star HD84937 is constrained to 
Dq < 4; this implies that at least one star on the Spite plateau has not depleted its 
primordial ^Li by more than a factor 4, even by extreme dilution. This constraint 
is in fact stronger when one takes all constraints into account; indeed, no current 
stellar model is able to reproduce the abundances of ^Li as a function of metallicity 
and effective temperature, and yet allow Dj > 2, while keeping Dq < 4. Therefore, 
'^Li should not be depleted by more than a factor ~ 2 on the Spite plateau. If direct 
nuclear burning were the depletion mechanism, then ^Li would be depleted by less 
than 2%. Moreover, all three ^Li observations are in excellent agreement with all 
standard expectations: Big-Bang nucleosynthesis with 2 < r/io < 6.5, and standard 
stellar isochrones of ^Li survival in metal-poor stars. We also discuss possible devi- 
ations from our assumption of a scaling of ^Li/H with due to q — q creation 
of ^Li in various sites. 
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1 Introduction 



The light nuchde ®Li is not significantly produced in standard Big-Bang nucleosynthesis, 
(^Li/H)p ~ 10~^^ (Thomas et al., 1993; see below for a new evaluation), and is expected 
to be produced over the hfetime of the Galaxy in galactic cosmic rays (GCR) p,Q;-C,N,0 
spallation, as well as a — a fusion reactions (Epstein, Arnett & Schramm, 1974; Reeves, 
1994, for a review). Its high fragility to stellar processing makes it a powerful tool to 
constrain Big-Bang nucleosynthesis (Epstein, Arnett & Schramm, 1974), as was pointed 
out by Brown & Schramm (1988). It was argued that the survival fractions ge, gr, of ^Li 
and ^Li, defined as the ratios of the observed to the proto-stellar abundances {i.e. =1/D, 
where D is the depletion factor), should be related by g^ a; g^, where /9 ~ 60 — 76 is 
the ratio of the destruction rate of ^Li to that of ^Li, in the case of nuclear destruction. 
Since g-j = 0.99 already corresponds to a factor ~ 2 destruction for ^Li, and since the 
expected proto-stcUar abundance of '^Li in metal-poor stars is expected to be at the limit 
of feasible detection (see below), the detection of ^Li in a metal-poor halo star could 
then be turned into a strong qualitative argument for the absence of significant (nuclear) 
depletion of ^Li. In that case, the uniformity of the ^Li abundances observed in metal- 
poor dwarfs on the so-called Spite plateau (Spite & Spite, 1982) should indeed reflect the 
primordial abundance of ^Li. Such an abundance, (^Li/H)p ~ 1.6 x 10~^°, would be in 
excellent agreement with standard Big-Bang nucleosynthesis and the inferred primordial 
abundances of D and ^He (Copi et al., 1995, and references therein). It should be noted 
that models where the surface abundances are depleted by mixing, and thus dilution with 
regions devoid of either ^Li of ^Li, are limited directly by the depletion factors. 

It is extremely difficult to detect ®Li in population II (PopII) stars, due to the weak iso- 
topic separation of the ^Li-^Li A6708A hues, as compared to the width of the ^Li line, and 
to the low expected isotopic ratio, ^Li/^Li^ 0.1, notwithstanding any possible destruction 
of ®Li. A first detection was nonetheless reported by Smith et al. (1993) in HD84937, 
^Li/^"'"''Li=0.05±0.02. This result was confirmed independently by Hobbs & Thorburn 
(1994), who reported 6Li/^Li=0.07±0.03 in the same star, as weU as 6Li/^Li=0.05±0.02 
in HD201891; another marginal detection was reported in HD160617 by Nissen (1995), 
6Li/6+7Li=0.017±0.012. 

Smith et al. (1993) were quick to point out the significance of their first measurement 
in HD84937. Using the observed abundances of boron and beryllium at low metallici- 
ties, and the {Li/B)GCR production ratio expected from GCR spallation reactions in a 
PopII environment (Steigman & Walker, 1992), these authors showed that the isotopic 
ratio ^Li/^Li was in excellent agreement with: standard Big-Bang nucleosynthesis, i.e. 
(^Li/H)p ~ 10^^°, no depletion of ^Li, possibly little destruction of ®Li, and no signifi- 
cant primordial production of ®Li. Moreover, these authors showed that in the case of 
rotational stellar models, which predict large destruction factors for '^Li, g-j ~ 0.1, and 
for ^Li, gQ ^ 0.04, the measured isotopic ratio in HD84937 imphes either one of the fol- 
lowing: (i) a massive production of lithium isotopes in a — a fusion reactions in a very 
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early phase of the Galaxy; (ii) a. very high primordial abundance of ^Li and ^Li, in a non- 
standard Big-Bang nucleosynthesis scenario; or (in) in- situ synthesis of ^Li in HD84937, 
for instance in stellar flares. We examine these different cases in Sec. 3. Steigman ct al. 
(1993) reached very similar conclusions using the observed ^Be abundance in HD84937: 
they found agreement with standard Big-Bang nucleosynthesis, with no depletion of ^Li 
and ^Be, and possibly moderate depletion of ^Li, qq ^ 0.2, as predicted by non-rotating 
standard stellar models. They also argued that these ^Li observations posed a severe 
challenge to the so-called Yale rotational stellar models. 

A potential problem with these two approaches is that they rely entirely on the knowl- 
edge of the (Li/Be)GCR production ratio at low metallicities, and it is now accepted that 
this ratio is in fact very uncertain (Fields et al., 1994, 1995). Indeed, the naive extrapola- 
tion at low metallicities of standard galactic cosmic ray (GCR) spallation is in apparent 
conflict with the observed abundances of ^Be and B in PoplI stars, since it predicts a slope 
~2 for the correlations of Log(Be/H) and Log(B/H) vs [Fe/H], whereas a slope ~ 1 ± 0.1 
is observed {e.g. Rebolo et al., 1988, Gilmore et al., 1992a, b, Ryan et al., 1992, Boes- 
gaard & King, 1993, and references, for Be; Duncan et al., 1996a, b and references, for B). 
In order to explain these primary behaviors of ^Be and B, different scenarios have been 
suggested: some involve modifying the GCR scenario (Duncan et al., 1992; Prantzos et 
al., 1993; Fcltzing & Gustafsson, 1994), some others modify the chemical evolution in the 
halo phase (Ryan et al., 1992; Fields et al., 1993, Tayler, 1995), while some others invoke 
different spallation processes (Casse et al., 1995; Vangioni-Flam et al., 1995). Moreover, 
whereas Steigman & Walker (1992) have shown that the (Li/Be)G'ci? production ratio 
depends in a crucial way on the metallicity of the medium, Fields ct al. (1994, 1995) have 
shown that it also depends quite strongly on the physical parameters of the GCR flux 
{e.g. spectral exponent and path length). Overall, it appears that this production ratio 
can be very different in the above models of early ^Be and B evolution, taking values in 
the range ~ 10 - 1000. 

In the present paper, we offer a simple prescription that allows us to circumvent the 
model dependent uncertainties, hence to derive empirical proto-stellar ^Li abundances. 
We simply assume that ^Li scales as oxygen at all metallicities. Such a behavior is indeed 
observed for ^Be, and any (p,Q;-C,N,0) spallation process which creates ^Be should also 
produce ^Li, in the same way, and vice- versa (sec also below). We justify this prescription 
further in Sec. 2. We then construct a predicted evolutionary curve for ^Li, and normalize 
it to the ^Li meteoritic abundance. It is then straightforward to obtain the possible 
depletion factors for ^Li in HD84937, HD201891, and HD160617 (Sec.2). We also discuss 
alternatives to this model (Sec.3): these can only arise from a — a fusion, since this is the 
only channel whereby ^Li and not ^Be can be produced. 

2 Standard expectations 
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(1) 



where i runs over p, a, and Z runs over C,N,0; y denotes the abundance with respect 
to hydrogen; the average is performed over energy and source and target composition; 
stopping factors due to energy losses and/or escape are implicitly included as form factors 
on the cross-sections. As has been extensively discussed in the literature (Vangioni-Flam 
et al., 1990; Fields et al., 1994), the standard GCR yields are proportional to the time 
derivative of the oxygen abundance squared i/q, giving rise to the slope 2, because the 
GCR flux (pi is proportional to the supernova rate, the yz abundances are proportional to 
the integrated supernova rate, and the flux 0^ is proportional to the product of these two 
factors. As discussed in Fields et al. (1994) and Lemoine et al. (1996), one can escape 
this slope 2 in at least three ways: (ijhj adjusting the Z dependence of the flux 0j (Fields 
et al., 1993); (ii) by keeping yz constant (Gilmore et al., 1992b; Feltzing & Gustafsson, 
1994; Tayler, 1995); and (Hi) by adjusting the (pz dependence on the supernova rate 
(Duncan et al., 1992; Vangioni-Flam et al., 1995; Casse et al., 1995). The inclusion of 
Tayler (1995) in category (ii) is not obvious, due to the peculiarity of this model; it is 
discussed in Sec. 3. Note that Feltzing & Gustafsson (1994) have shown that spallation 
in the vicinity of supernovae, as advocated by Gilmore et al. (1992b), is probably not 
energetically viable, hence we will refer to such models only as an example of category (ii) 
above. As well, we note that Lambert (1995) proposed to explain the halo abundances of 
^Li, Be and B by direct production of these elements in the observed stars from cosmic 
rays impinging on the stellar atmospheres. He showed that this scenario was not viable 
for Be and B, and barely able to reproduce the ^Li abundances. 

It is then straightforward to derive the scaling of ^Li with metallicity, rewriting Eq. (|l]) 



In cases (ii) and ( Hi ), the yields of ^Li scale as the supernova rate. Only in the first scenario 
does the behavior of ^Li deviate from that of a primary element. Indeed, appropriately 
adjusting the flux can make the ^Be correlation go from quadratic to linear, in which case 
the ^Li correlation goes from linear to ~constant through the a — a fusion channel. This 
is easily seen by parametrizing the metallicity as a power-law in time, ype ^" (note 
that ype is the iron abundance not in logarithmic units); tuning the flux so as to obtain a 
slope 1 for Be and B then leads to: log(^Li/H)oc log (constant + [Fe/H]). However, as we 
discuss in Sec. 3, it is not obvious why and how such a tuning of the flux could take place 
in the halo phase. 

This brief discussion shows that, under very general conditions, ^Li should behave as 
a primary element in the halo phase. However, our prescription further requires that this 



for "Li: 



dyu 



{(paCT^aVa) + {(piCrl^^z) + {(pZCT^^iVi) ■ 



(2) 



dt 
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scaling holds equally well up to solar metallicities, i.e. that the transition from the halo 
to the disk phase proceeds without any change in the (^Li/0) production ratio. Here, 
we distinguish two different cases: (i) where GCR are responsible for the halo evolution 
of LiBeB, through some modification of their parameters or some modifications of the 
galactic evolution; and (ii) where another process than GCR accounts for the evolution of 
LiBeB in the halo phase. In the first case (i), we note that, for GCR spallation and fusion, 
at a given (low) metallicity, although (^Li/0) depends more strongly than (^Be/0) on 
the very low energy part ~ 10 MeV/n of the GCR spectrum, due to the low energy a — a 
cross sections, it depends less strongly on the spectral slope and the path length, that 
affect the higher energy part of the spectrum. We note as well, that, within observational 
accuracy, the (^Be/0) ratio keeps constant from the halo to the disk phase (see Fig.H). It 
is therefore reasonable to assume that, in case (i), the (^Li/0) production ratio has not 
changed significantly between the halo and the disk phase. In the second case (ii), we 
can reasonably assume that the process responsible for the evolution of Be and B in the 
halo phase, proceeds on into the disk phase, and that the production ratio (^Li/0) ratio 
does not change. This is indeed what Fig.|^ suggests. One knows as well that standard 
GCR spallation/fusion creates ^Li as a primary element in both halo and disk phase. 
When mixing these two processes of creation of ^Li, one mixes two primary processes, 
and one obtains a primary process. In this case, therefore, there is no change of slope or 
of production ratio (^Li/0). 

Finally, the constancy of the ratio (^Li/0) ratio is shown to hold for the models 
referred to in (ii) and (Hi), except Tayler (1995); in this latter scenario, the (^Be/O) ratio 
in the halo phase is arbitrarily tuned to that in the disk phase, but it is not possible to 
make definite predictions for the (Li/0) ratio. We discuss this model further in Sec. 3. In 
any case, we allow, in assigning error bars for our predicted evolutionary curve, for some 
possible (slight) mismatch in the (^Li/0) production ratio between the halo and the disk 
phase. 

2.2 ^Li as a primary element 

In order to construct the evolutionary curve of ^Li, we use a simple model of chemical 
evolution. This model is a closed box, one zone model based on the Tinsley (1980) 
formalism, able to reproduce the trends in C, N, O, [C/Fe]~[N/Fe]~ 0, [0/Fe]~ 0.5 in 
the halo phase (Wheeler et al., 1989), together with the age- metallicity relation, and the 
solar abundances of C, N, O, and Fe at solar metallicity. This model incorporates a GCR 
code for production of LiBeB. We used this model in order to take into account the slope 
changing effect on the ^Li-Fe relation in the disk phase, due to the introduction of Type 
la supernovae, i.e. in order to follow [0/Fe]. It is sufficient for our purpose to use the 
standard GCR calculations of spallation/fusion production of ^Li, with a flux proportional 
to the supernova rate, as these calculations yield the desired primary behavior. More 
complex issues of galactic evolution, such as supernovae ejecting metals to form a hot 
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intergalactic gas (Mathews & Schramm, 1993) are irrelevant here. 

We note that the upper hmit on the primordial abundance of ^Li has recently been 
revised upwards by two orders of magnitude, through new measurements and limits on 
the cross-section d(a,7)^Li at low energies (Cecil et al., 1996). We included this upper 
limit in the predicted curve, viz. (^Li/H)p < 1.2 x 10~^^, although it has no effect on our 
calculations, since the metallicities where ^Li has been detected are too "high" for this 
primordial component to be relevant. Nonetheless, as we show in Fig.|l], this could be 
of relevance for metallicities [Fe/H]^ —3, depending of course on the actual primordial 
abundance of ^Li. This will be discussed in detail in Nollett et al. (1997). 

We thus construct the ^Li evolutionary curve, and normalize it to the meteoritic 
abundance of ^Li, as shown in Fig.^. Data points for ^Li were taken from Smith et 
al. (1993), Hobbs & Thorburn (1994), Nissen (1995), and Nissen et al. (1995), and 
are reproduced in Table |l]. We did not consider in these data upper limits that were 
previously obtained by Andersen et al. (1984), Maurice et al. (1984), Hobbs (1985), and 
Pilachowski et al. (1989); although these observations certainly represented a challenge 
at that time, the upper limits obtained, ^Li/^Li< 0.1 in all cases, are unfortunately not 
stringent enough for our purpose. Data points for ''Li and ^Be come from a compilation 
of Spite & Spite (1982), Boesgaard & Tripicco (1986), Hobbs & Duncan (1987), Rebolo et 
al. (1988), Pilachowski et al. (1993), Gilmore et al. (1992a, 1992b), Ryan et al. (1992), 
Boesgaard & King (1993), Smith et al. (1993), Hobbs & Thorburn (1994), Thorburn 

(1994) , Duncan et al. (1996a, b), Molaro et al. (1995), Nissen et al. (1995), Rebolo et al. 

(1995) , Spite et al. (1996), and Anders & Grevesse (1989). 

The ^Li abundances, [Fe/H], and T^ff for the ^Li stars were taken from Smith et al. 
(1993), Hobbs & Thorburn (1994), Thorburn (1994), Nissen (1995), Nissen et al. (1995), 
Primas (1995), Spite et al. (1996), and Ryan et al. (1996); see Table We adopted, as 
error on the temperatures, 150K statistical, and lOOK systematic, at 95% c.l., following 
Spite et al. (1996). Note that, contrary to elemental abundances, the isotopic ratio is 
only weakly sensitive to the [Fe/H] and T^/f adopted; for instance, the ^Li abundance of 
HD84937 as derived by Thorburn (1994) is 0.24dex higher than that derived by Smith 
et al. (1993), although the isotopic ratios are entirely consistent. For this star, both 
®Li/®"''''Li evaluations were combined in a weighted least squares procedure. 

We adopted a ±0.16dex 2a statistical gaussian uncertainty for lithium and iron abun- 
dances, which is a typical estimate (Thorburn, 1994; Spite et al., 1996). We adopted a 
±0.2dex 95% c.l. systematic uncertainty for lithium and iron abundances, a typical esti- 
mate as well (Thorburn, 1994; Spite, 1995; Spite et al., 1996). In Fig.^, we combined these 
systematics on the ''Li and ^Li abundances quadratically with the 2a statistical error bars, 
in order to produce an estimate of the total 95% c.l. on these data points. In the following 
analysis, however, we allow for different statistics. Namely, we take the abundance sys- 
tematics on N[Li] and [Fe/H] to follow top-hat probablity laws, and the ^Li/^"^^Li ratios 
to follow gaussian laws. We adopted a ±0.3dex 95% c.l. uncertainty on our evolutionary 
curve, which we can attribute to the normalization to the meteoritic abundance of ®Li 
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as well as to slope changing effects of chemical evolution at the halo-disk transition. In- 
deed, although the meteoritic abundance of ^Li, as given by Anders & Grevesse (1989), 
is endowed with very small error bars, ±0.04dex, we prefer to remain conservative and 
allow for a possible la uncertainty ~ ±0.2dex in promoting this meteoritic abundance to 
a 'cosmic' abundance at [Fe/H]=0. We show, in Fig.^, how well this evolutionary curve, 
and its accompanying error bars, bracket the Be observations in halo stars. We always 
refer to ^Be as a baseline, and not to B, since the available boron abundances have not 
yet been fully corrected for NLTE effects. Let us mention, however, that a clear primary 
correlation for boron vs iron in the halo phase, with a slope ~ 1.0 ± 0.05, is inferred from 
these latter abundances (Duncan et al., 1996a,b). 

We now discuss the three detections of ^Li in HD84937, HD201891, and HD160617. 
The upper limits located above the evolutionary curve simply represent non-detection of 
^Li, while those below the curve are interpreted as non-detections due to destruction of 
^Li in these stars. We will show in later discussion that this is compatible with their Tejj. 
It is straightforward to evaluate the survival fraction of ^Li in the three stars through 
statistical bootstrapping, using the estimator: 

In using this estimator, it is implicitly assumed that no production of ^Li took place 
in the star, i.e. we neglect any flare production (Deliyannis & Malaney, 1995); we will 
return to this point in Sec. 3. In this bootstrapping, we draw at random, according to their 
statistics, the 'observed' and the 'predicted' abundances. We draw the predicted proto- 
stellar abundances for the three stars from the same evolutionary curve, so as to ensure 
that the slope of the ^Li correlation vs [Fe/H] is 1 in the halo phase. We do not correlate 
the systematics of the observed abundances, since these stars have different metallicities 
and different Tg/j. We took care to include the metallicity statistics of the stars in the 
bootstrapping, taking into account the correlation of (®Li/H)pred with metallicity. The 
resulting qq histogram is shown in Fig.^ for a sample size of 10^. 

We thus obtain the following averages and 95% c.l. limits on (yfg, in the standard case: 

HD84937 = 1-92 0.57 < < 7.24 
HD201891 = 0.14 0.02 < < 0.54 
HD160617 gQ = 0.38 0.00 < < 1-67. 

The lower limit on g^ for IID160617 is 0.00: it could have been expected on the basis 
that ^Li was not detected at the 2a level in this star. We checked that the figures above 
were not sensitive to the type of statistics adopted for the systematics and the theoretical 
modeling, within ~10%. 

With our present assumption that flare production of lithium did not occur in these 
stars, > 1 is unphysical. One usual remedy is to use Bayesian statistical inference. 
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and to truncate the probability with (for instance) a top hat function 9 {g^) 6 (1 — g^), 
to enforce < ge < 1, and re-normahze the cjq distributions; in that case, the resulting 
95% c.l. limits would become (0.26,0.99) for HD84937, (0.03,0.48) for HD201891, and 
(0.03,0.91) for HD160617, and the averages: 0.71 for HD84937, 0.14 for HD201891, 0.30 
for HD160617. We note, however, that these numbers are subject to the (subjective) 
choice of constraints applied on g^, as is the case in any Bayesian inference procedure. 
The use of a top-hat uniform law is justified as a zero-th order choice when one does not 
have any more reliable guess to offer. This is the case here, since so-called 'standard' 
isochrones do not agree with each other as to the degree of destruction of ^Li in these 
stars (see Chaboyer, 1994, on the one hand, and Deliyannis et al., 1990, on the other 
hand) . 

Another remedy is to invoke potential uncontrolled systematics, which, in the present 
case, should be associated with the estimates of the ^Li/^+^Li ratio. Indeed, as pointed 
out by Nissen (1995), convective motions in the stellar atmosphere result in an asymmetry 
of the profile, that might mimic ^Li absorption in the Lil line. The choice of the type 
of absorption profile also plays an important role in the estimate of the ^Li/^"^''Li ratio. 
Although these effects were taken into account by Nissen (1995), Nissen et al. (1995), 
using the nearby Fel and Cal lines as standards, the choice of these lines is dictated by 
convenience alone, and some mimicking of ^Li is still possible. Such systematics could 
be ascertained using a greater sample of lines as standards; obviously, higher spectral 
resolution at high signal-to-noise ratio would allow a substantial increase in the accuracy 
of these estimates. A zero-th order estimate of such systematics is given by the difference 
of the observed to the predicted abundance of ^Li in HD84937, yielding a reasonable offset 
of -|-0.3dex. This star is indeed the hottest of the three stars, and, for T^jf ~ 6300, possibly 
no depletion of ^Li is expected in stellar models (Brown & Schramm, 1988; Deliyannis et 
al., 1990, Chaboyer, 1994), depending on the (unknown) evolutionary status of HD84937, 
whose photometry puts it at the turn-off. Note that such systematics would not be in 
contradiction with the other two stars, that are main-sequence dwarfs with Te// ~ 5800, 
for which moderate to strong ^Li depletion is expected. Thus, if we shift, by hand, the 
observed data point, to place it right on top of the predicted curve, the 95% c.l. limits 
on in HD84937 become (0.30,3.16). It is comforting to note that both remedies used 
above, to enforce < 1, lead to a very similar 95% c.l. lower limit on gg. 

Finally, in order to check the consistency of the derived depletion factors with stellar 
isochrones for ^Li survival, we plot in Fig.^ the observed data points and upper limits vs 
Teff. Unlike ''Li, however, for which the measured abundances are at the same level on 
the Spite plateau, the ^Li data points need to be corrected for their trend with metallicity; 
in Fig.^ they were brought back to the metallicity of HD84937. It is obvious from Fig.^ 
that the locus of these data points is consistent with the expected shape of isochrones, as 
shown in Chaboyer (1994) for instance, and, more specifically, that none of the upper limits 
derived contradict the three detections. Moreover, we did not reproduce the isochrones 
of Chaboyer on Fig.^ since most of these stars are located near the turn-off, and the 
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isochrone seems to depend strongly on the evolutionary status of the star (Chaboyer, 
1994). As well, note that to change the opacities used in the code is sufficient to modify 
the ^Li depletion factor from to ~ 3 (in the subgiant case), which means, conversely, 
that we cannot rely too much on these models. In any case, the general agreement is very 
satisfying. 

To conclude this section, we obtain a 95%c.l. lower hmit on the survival fraction 
in HD84937, > 0.26 (equivalently Dq < A for the depletion factor), assuming some 
uncontrolled systcmatics in assessing the line profile of the ^'^Li A6708A line, or using 
Bayesian inference to bring back the measured average in the physical region < 1. 
It is well known that ^Li is much more fragile to nuclear burning than ^Li {e.g. Brown 
& Schramm, 1988, and references therein). As a consequence, the preservation region of 
^Li is, at any time, much shallower than that of ^Li. As a consequence, ^Li is also more 
sensitive than ^Li to depletion mechanisms such as dilution processes, or even mass-loss 
(Vauclair & Charbonnel, 1995). Therefore, it seems reliable to use the lower limit on gQ 
as an extreme lower limit on the ^Li survival fraction g'^. Our evaluations thus show that 
at least one star on the Spite plateau has depleted its ^Li by less than a factor 4, at 95% 
C.I.. Using the observed abundance of ^Li in HD84937, Log(^Li/H)=-9.77 ± 0.16 ± 0.2 
(95% c.l. estimates for both statistical and systcmatics), we derive an upper limit on the 
primordial abundances of ^Li, Log(^Li/H)p < —8.81 (95% c.l.), where we added linearly 
the errors and the upper limit on the depletion factor. Although this upper limit is not 
too stringent in itself, we note that, in order to saturate the upper bound, one would 
have to deplete both ^Li and ^Li by the same amount (a factor 4), and yet reproduce all 
the observed features of the Spite plateau, namely the uniformity in abundances and the 
isochrones. The above constraint is therefore stronger, i.e. it amounts to D-j ^ 2, since 
there is, at our knowledge, no stellar model able to reproduce the observed trends, while 
having > 2 and Dq < 4. We note, for instance, that, if the maximal depletion of ^Li 
were due to pure nuclear burning without dilution, then ^Li would be depleted by less 
than 2% from its primordial value (Brown & Schramm, 1988). As well, the rotational 
mixing models (Pinseonneault et al., 1992) satisfy the latter constraint on the features of 
the Spite plateau with a primordial abundance Log(^Li/H)p ~ —9, but cannot do without 
a depletion factor Dq ^ 25 (Dehyannis & Malaney, 1995). 

3 Alternative scenarios 

As mentioned earlier, the only way ^Li can avoid a primary behavior for the whole range 
of metallicitics —4 <[Fc/H]< 0, comes through the a — a fusion channel. To discuss these 
'non-standard' cases, we are compelled to restrain ourselves to specific models: arbitrary 
tuning of the cosmic ray flux, referred to as (i) in Sec. 1.1, spallation in globular clusters 
(Tayler, 1995), or stellar flare production of ^Li (Deliyannis & Malaney, 1995). 



9 



3.1 Early bright galactic phase 

As discussed earlier, different models have been proposed to tune the cosmic ray flux so 
as to reproduce the observed abundances of ^Be and B in the galactic halo. Ryan et al. 

(1992) suggested an outflow model with an enhanced supernova rate in the halo phase; 
since the injection flux is traditionally assumed to be proportional to the supernova rate, 
it was believed that it would allow the production of higher Be and B abundances. This 
is incorrect in that increasing the supernova rate leads to an increase in the production 
of metals, so that a higher ^Be or B abundance would be associated with a higher metal 
abundance, and no difference with the standard model would be noted in an elemental 
graph Log(^Be/H) vs [Fe/H]. The higher ^Be and B abundances were rather obtained in 
an ad-hoc fashion by normalizing the yields to the observed ^Be abundance at [Fe/H]=— 1; 
in particular, a slope 2 is obtained. The role of the outflow was to evacuate the gas to 
prevent a high supernova rate in the disk phase, so as to prevent an overproduction of 
^Be and B at solar metallicities (see Fig.3 in Prantzos et al., 1993 for an illustration of 
these effects). Therefore, since further observations of Pop 11 ^Be and B abundances seem 
to conflrm a slope of 1 on the whole range of metallicities —4 <[Fe/H]<0, we need not 
consider this model further. This remark also applies to the models of Prantzos et al. 

(1993) , where the path length was tuned in the halo phase so as to reproduce the ^Be and 
B abundances; although this model was nicely motivated, namely the cosmic rays should 
be overconflned in the halo phase, it led to a slope ~ 1.75 for ^Be and B, and, in any case, 
to a slope ~ 1.0 for Li. 

Finally, in Fields et al. (1993) it was shown that to reproduce a slope of 1 for ^Be and 
B, the cosmic ray flux should vary approximately inversely with time. These authors did 
not offer an explanation for such a behavior: they were rather interested in constraining 
the duration and luminosity of an early bright galactic phase under the assumption that 
the cosmic ray flux would effectively vary as 1/t. It was shown that the bulk of the 7-ray 
background could be associated with such an enhanced cosmic ray flux in the flrst few 
bilhon years. Nevertheless, according to the above discussion, Ryan et al. (1992) have 
shown that an increased supernova rate would not lead to a slope 1, and Prantzos et 
al. (1993) have shown that an overconflnement would not do so either. Hence, the only 
remaining explanation for a GCR flux varying as the inverse power of time would be to 
have the injection flux evolve as a high power of the supernova rate, which does not have 
a physical basis at the present time. Indeed, with (0/H)oc i", n ^ 1, one would require 
(f)inj oc (dNsN / dt)^^^~^ . Finally, we note that Lemoine et al. (1996) have shown that 
the GCR propagated flux should evolve as a very low power < 1 of the supernova rate, 
throughout the disk phase, in order not to overproduce the LiBeB abundances at solar 
metallicity. Thus, the higher power could only apply in the early pre-disk phase. 
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3.2 Spallation in globular clusters 

Tayler (1995) has offered an original explanation to the primary behaviors of ^Be and B in 
the halo phase, in suggesting that the bulk spallation would take place in globular cluster 
sized objects before the gas would be swept out of the cluster and mixed with the ISM 
gas. The primary behavior arises from the assumption that every globular cluster ejects 
the same amount of gas, and that the spallation is the same in every cluster: each time 
metals (Z) are ejected to enrich the ISM, LiBeB elements are also ejected, with a constant 
(LiBeB/Z)ejec ratio, that is independent of the Z content of the ISM. This ratio should not 
in fact be constant, as these globular clusters should not eject exactly the same amount 
of metals, and the spallation is not expected to be exactly the same in every object, so 
that, in the end, one would expect to observe a linear relation Log(^Be/H) vs [Fe/H] with 
some scatter. There has indeed been a recent claim for a scatter around the linear trend 
(Primas, 1995), which gives credit to this scenario. 

This model is nonetheless subject to relatively strong requirements. It assumes that all 
globular clusters eject nearly the same amount of metals, which, in a first approximation, 
would mean that they all have very similar masses. Moreover, it is assumed that the 
spallation is much more efficient in the globular clusters than in the ISM. Finally, the 
(^Be/0) production ratio has to be arbitrarily tuned to that observed in the disk phase, 
since the linear behavior of ^Be vs [0/H] is observed on the whole range —3 ^[O/H]^0 
(Fig.|]). On this basis, it is unfortunately difficult to say whether or not the (Li/0) 
production ratio is the same in these globular clusters and in the disk phase, since the 
(Li/Be) ratio depends very strongly on the spectral slope and the confinement time of the 
accelerated projectiles (Fields et al., 1993, 1995). Hence, no reliable evolutionary curve 
for ^Li can be constructed in this case. 

3.3 Flare production 

Deliyannis & Malaney (1995) have studied the possibility of generating the observed ®Li 
of IID84937 via flare production. They argue that the energetics can, in principle, be 
fulfilled. However, we disagree with their estimate of the energetics required, and we 
show below that such a model is strongly disfavored for producing all the observed ^Li. 
These authors also propose two observational tests to possibly discriminate in situ flare 
produced ®Li from proto-stellar ^Li: (i) the ^Li/^Be and B/Be ratios in HD84937; and 
(a) the curve of ^Li/^Li ratios the effective temperature in a sample of metal- deficient 
dwarfs. We disagree with these authors on the feasibility and the reliability of these tests, 
as we discuss now. 

The first test of Deliyannis & Malaney (1995) to distinguish in situ flare produced 
®Li from proto-stellar ^Li relies on the difference between the ^Li/^Be and B/Be ratios 
resulting from flare production and those resulting from cosmic ray spallation. However, 
it is straightforward to show that, even if all ^Li observed in HD84937 were produced by 
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flares, the amount of Be and B produced by these flares would be neghgible with respect 
to the amount of Be and B observed in this star. Using the flare production ratios of 
Walker et al. (1985), ^Li/^Be~ 40 — 100 for a solar metallicity star, and rescaling it to 
a metallicity [Fe/H]=-2.1, this ratio should become ^Li/^Be- 10^ - 10^ in HD84937, 
since ^Li is still produced through a — a fusion, but the p-C,N,0 channel is cut off. This 
ratio has to be compared with the observed ratio ^Li/^Be~ 66^4° (95% c.L): clearly, if 
all observed ^Li has been produced in situ by flares, then only a tiny (and negligible) 
fraction of the observed ^Be comes from these flares. This conclusions holds true for B 
as well, and therefore, the ^Li/Be and B/Be ratios cannot be used to trace in situ flare 
production. 

The Dehyannis & Malaney second observational test rests on the dependence of the 
observed ^Li/^Li ratio with the effective temperature in subgiants. Deliyannis & Malaney 
(1995) have shown that this ratio should remain constant if the observed ^Li is of proto- 
stellar origin, in the range of temperatures 6400> Tg// > 5800, as neither of both isotopes 
is then depleted. If ^Li has been created in flares, the isotopic ratio should decrease with 
decreasing temperature, in the same range of temperatures, due to increased dilution 
(i.e. larger convection zone) past the turn-off. Indeed, the essential difference between 
proto-stellar and flare produced ^Li is that the proto-stellar ^Li preservation region is con- 
siderably larger than the convection zone post turn-off, whereas the only region containing 
flare produced ^Li is the convection zone itself. This behavior is, however, very dependent 
on the evolution of the convection zone, which is difficult to model during the ascent on 
the subgiant branch. Moreover, as pointed out before, it seems that the survival of ^Li in 
these stars depends strongly on the opacities adopted. Overall, it seems that, even if a 
large database of ''Li/^Li ratios in metal-poor stars could be collected, a possibility that 
is already very optimistic, this would not readily allow one to discriminate between flare 
produced and proto-stellar ^Li. 

We now turn to the energetics required for producing ^Li in flares in HD84937. Deliyan- 
nis & Malaney (1995) argue that the production efficiency is 10~^-1 atom/erg, so that 10 
flares per year of strength 10"^^ ergs, which corresponds to the large solar proton flares, 
for 1 atom/erg, during IGyr, would reproduce the observed ^Li. The point is that large 
proton flares, in the Sun, are associated with steep power law spectra (Ramaty et al., 
1995), and the efficiency for these spectra is ^ 0.007 atom/erg (with a low-energy cut-off 
imposed at 1 MeV) for a solar composition (R. Ramaty, 1996, private communication). 
We calculate that for a composition rescaled to the metallicity of HD84937, this estimate 
goes down to ~ 0.003 atom/erg, assuming a 1 MeV cutoff. The Deliyannis & Malaney 
estimate of 1 atom/erg is therefore too optimistic by a factor ~ 300, and instead of 10 
flares per year, ~ 3000 flares per year of a strength comparable to that of the large solar 
flares would be required, an enormous number indeed. 

Here we investigate in greater detail the production of lithium in flares. We use only 
the a — a fusion channel in our calculations, since with its threshold and peak resonances 
around ~ 10 MeV/n, and the low metal content of the star, it should largely dominate 
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over the p,a-C,N,0 channel. We assume an abundance na/rip = 0.1 in the accelerated 
and target particle abundances, although for large solar flares, the accelerated abundance 
could be lower by a factor 10 (Murphy et al., 1991). We also assume that, through 
downward convection, all the ^Li produced in flares would be mixed in the convection 
zone of mass Mc. In a thick target model, the rate of ^Li production can be written 
(Ramaty & Murphy, 1987): 

where Q (^Li) is in atom.s"^, {dE/dx)^ represents the energy losses of a nuclei in the in- 
teraction medium, in MeV.g~^.cm^, and Na [E) is the injection spectrum of a particles in 
(MeV/n)-^s-^ This latter can be parametrized by a power law, a decreasing exponential 
in rigidity, or a Bessel function K2 of the second kind, although we note that this latter 
reproduces more adequately the observational data on solar flares. It is also customary 
to define the deposited power: 

r+00 

W = J2A EN,{E)dE, (5) 
Jo 

where Ai denotes the atomic mass of nuclei de i. The production efficiency, referred to 
above, is then defined as Q (^Li) /W, in atom.erg"^. Here, however, we calculate the 
rate of production of ^Li nuclei, using Eq.§, normalizing the injection spectrum to the 
observed rate of average irradiation of the solar atmosphere, i.e. Np (> 30MeV) ~ 3.10^^ 
protons. s~^ (Ramaty & Simnett, 1991). We find that the production rate does not depend 
strongly (within a factor 2) on the type of spectrum chosen, or on the low-energy cut- 
off adopted for power law spectra. For power law spectra, we adopt a steep spectrum 
(index ~ 4), with a low-energy cut-off of 1 MeV/n. Apart from preventing divergence, 
this cut-off is to simulate the presence of energy losses that flatten the spectrum at low 
energies. Due to the normalization adopted above, and the fact that the threshold of the 
a + a ^ ^Li cross-section is ~ 10 MeV/n, our choice of cut-off maximizes the production 
rate. The production efficiency is strongly dependent on the type of spectrum, varying 
between ~ 10~^ and ~ 10~^ atom.erg"^, in agreement with previous calculations (Canal 
et al., 1975, 1980). The production rate obtained is: Q (^Li) ~ 2 - 4 x 10^^ atom.s'^ 
which, after mixing in the convective zone, on a production timescale r = 1 Gyr, yields: 

- 2-4 X 10" 

As compared to the observed ^Li abundance in HD84937, (^Li/H)~ 9. x 10^^^, this 
estimate is negligible. It has to be noted, however, that the mass Mc of the convection 
zone is largely unknown, mostly because the evolutionary status of HD84937 is itself 
unknown. Deliyannis & Malaney (1995) point out that Mc should range from 3 x 1O~^M0 
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to 3 X IQ-^Mq, with an estimated IQ-^Mq, if HD84937 is a dwarf, and from 7 x IO'^Mq 
to 1 X 10~^Mq, with an estimated 3 x 10"^ if HD84937 is a subgiant. The mass of the 
convection zone also depends strongly on the stellar parameters, such as the mixing length, 
and, overall, the above estimates should be good within an order of magnitude. The 
estimate of a production timescale r ^ 1 Gyr stems from the dependence of the convection 
zone on age: since the convection deepens with age, on the main sequence, any spallated 
®Li is more likely to have been produced recently; Deliyannis & Malaney (1995) estimate 
1 Gyr to be a reasonable timescale, and we adhere to this value. This shows that, even 
in the most favorable case (for flare production), i.e. a convection zone Mc ~ 10~^ Mq, 
the ^Li abundance produced in flares remains more than an order of magnitude below 
the observed abundance. Recall as well that we assumed that the downward convection 
of flare produced ^Li was fully efficient, that we assumed a "maximal" Ua/rip abundance, 
and that we maximized the production rate in our choice of spectrum. 

Hence, we feel that flare production of ^Li cannot produce a significant contribution to 
the observed ^Li abundance, unless HD84937 is subject to tremendously powerful flares 
over more than 1 Gyr. 

4 Conclusion 

We have estimated the depletion factor of ®Li in the three metal-poor halo stars where it 
has been detected. In order to determine the proto-stellar ®Li abundance in these stars, 
we have constructed an "empirical" evolutionary curve for ®Li, deduced from the observed 
abundances of ^Be. Namely, we have assumed that ^Li/H scales as over the galactic 

lifetime, and normalized the resulting evolutionary curve to the meteoritic abundance of 
^Li. Our assumption is motivated by the recent observations of metal-poor stars, that 
revealed a slope ~ 1 for the correlation of Log(^Be/H) with Log(^^0/H). One may indeed 
recall that ^Li and ^Be share their origin in p,Q;-C,N,0 spallation reactions (apart from 
an extra a — a fusion channel for ^Li). 

We have taken care in assigning uncertainties to the so-modeled evolutionary curve 
of ^Li. Notably, we have included new cross-section data for the primordial production 
of ®Li. We form the ratio of the observed to the proto-stellar abundance to obtain the 
survival fraction ^(g, i.e. the inverse depletion factor; the statistics of this survival fraction, 
in each star, are obtained from bootstrapping. We thus obtain a 95%c.l. limit > 0.26 
(equivalently Dq < 4. for the depletion factor) in the hottest turn-off HD84937 {T^ff — 
6300K). We have argued that some systematics, arising from the inaccuracy of stellar line 
profile modeling (at the precision required), should be associated with these measurements 
of the ^Li/®"''''Li ratio, as there is always the possibility that some weak profile distorsion 
might mimic the ^Li absorption. The above upper limit has been obtained in two different 
ways: assuming such systematics on the ^'^Li line profile, and proceeding to Bayesian 
inference, to return the average ~ 2. (measured for HD84937) into the physical region 
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Since ^Li is much more fragile than ^Li to nuclear burning, and since the preservation 
region of ^Li is always shallower than that of ''Li, this upper limit on the depletion factor 

is also an extreme upper limit on Dj, the depletion factor for ''Li. At a minimum, this 
implies that at least one star on the Spite plateau, with abundance Log(^Li/H)=— 9.77 ± 
0.16 ± 0.2 (95% c.l. for each error bars, statistical and systematics) , has depleted its ^Li 
by no more than a factor 4. This, in turn, means that an extreme upper bound can be put 
on the primordial abundance of ^Li, Log(''Li/H)p < —8.81 (95% c.l.). We note that, in 
order to saturate this upper bound, ^Li and ®Li would have to have been depleted at the 
same level, by a factor 4. For instance, if depletion had occurred through direct nuclear 
burning, ^Li would have to be depleted by no more than 2% for ^Li to be depleted by 
a factor 4 (Brown & Schramm, 1988). As well, rotational mixing models (Pinseonneault 
et al., 1992), that claim to obtain a much higher depletion factor for ^Li than 'standard' 
models, and yet reproduce all the observed features of the Spite plateau, are ruled out: in 
these models, a primordial abundance Log(^Li/H)p ~ —9 is associated with a depletion 
factor Z^e ~ 25 — 40 (Deliyannis & Malaney, 1995). We know of no current stellar model 
able to reproduce the observed features of the Spite plateau, and yet allow > 2, while 
keeping < 4. The above constraint on the depletion factor of ^Li on the Spite plateau 
is therefore stronger, i.e. ^Li should not be depleted by more than a factor ~ 2. 

Finally, we showed that all three measured ^Li abundances, as well as the upper 
limits derived, are in excellent agreement with all standard expectations: standard Big- 
Bang nucleosynthesis with 2 < 7710 < 6.5 (Copi et al., 1995, and references), 'standard' 
stellar models with respect to ^Li and ^Li survival in metal-poor stars, and standard early 
galactic evolution of both lithium isotopes (primary yields in spallation-fusion processes). 
The graph of ^Li abundances vs T^jf obtained is also in good agreement with the expected 
shape for ^Li isochrones. 

It is obvious that future measurements of ^Li abundances at different metallicities 
would narrow down the uncertainty in the possible depletion factors for ^Li in these stars. 
It is equally obvious that such studies would greatly benefit stellar modeling physics and 
the determination of the primordial abundance of ^Li. We recommend measurements of 
the ®Li abundance at a metallicity [Fe/H]~ —1, since this the region where ®Li detec- 
tion is favored (see Fig.|l|). Moreover, such measurements, combined with those at lower 
metallicities, would allow the determination of the slope of the correlation Log(^Li/II) vs 
[Fe/H], which we assumed to be ~ 1. 

We discussed three alternatives to the above scenario. Such alternatives can only arise 
from a — a fusion creation of ^Li, since this is the only channel whereby ^Li and not ^Be 
is produced. A condition on such alternatives is that they have to respect the primary 
behavior of ^Be in the early Galaxy. 

In the first case, we considered a logarithmic slope for Log(^Li/H) vs [Fe/H], resulting 
from an arbitrary adjustment of the GCR flux to the inverse power of time. We argued 
that this scenario lacked any physical basis, since it would require the early GCR injection 
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flux to be tuned to a high power of the supernova rate, in order to reproduce the Be and 
B trends. 

In the second case, where the bulk of LiBeB spallation would take place in globular 
cluster sized objects (Tayler, 1995), no firm prediction as to the behavior of ^Li can 
be made at the halo-disk transition, hence no evolutionary curve can be constructed. 
This scenario remains as a possible loophole for obtaining higher depletion factors of 
^Li in HD84937. We note, in passing, that among the models proposed to account for 
the primary behavior of ^Be in the halo phase, only this one and that of Casse et al. 
(1995), Vangioni-Flam et al. (1995) seem credible at this time. In this latter, C,N, and O 
nuclei, freshly synthesized by supernovae exploding inside their parental molecular cloud, 
are accelerated and eventually decelerate or produce LiBeB nuclei through spallation on 
interstellar p and a. 

Finally, we showed that stellar flare production of ®Li cannot produce a significant 
contribution to the observed ®Li, contrary to the suggestion put forward by Deliyannis & 
Malaney (1995). 
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Table 1: Upper limits and measurements of the ^Li/^+^Li ratio in halo stars. Error bars 
are 95% c.l. estimates; when two error bars are quoted, the first corresponds to statistical 
errors, the second to systematics. References: (a) Ryan et al., 1996; (b) Thorburn, 1994; 
(c) Nissen, 1995; (d) Nissen et al., 1995; (e) Smith et al, 1993; (f) Hobbs & Thorburn, 
1994; (g) Primas, 1995; (h) Spite et al., 1995. Evolutionary Status: T-o stands for Turn- 
off, Dw for Dwarf, and Sg for Sub-giant. 



Star 




[Fe/H] 


N[Li] 


(•^Li/'^+^Li) 


Status 


Reference 


HD84937 


6280±150 ± 100 


-2.1±0.16±0.2 


2.23±0.16 ±0.2 


0.055±0.034 


T-o 


(a),(e),(f) 


HD201891 


5800±150 ± 100 


-1.4±0.16±0.2 


1.87±0.16±0.2 


0.048±0.040 


Dw 


N,(f) 


HD160617 


5800±150 ± 100 


-2.0±0.16 ±0.2 


2.11±0.16±0.2 


0.017±0.024 


Sg 


(c),(9),(h) 


BD 3°740 


6340±150 ± 100 


-2.9±0.16±0.2 


2.21±0.16 ±0.2 


< 0.08 


T-o 


(a),(f) 


BD 26"3578 


6150±150 ± 100 


-2.4±0.16±0.2 


2.12±0.16 ±0.2 


< 0.10 


T-o 


N,(f) 


HD19445 


5850±150 ± 100 


-2.1±0.16±0.2 


2.08±0.16±0.2 


< 0.03 


Dw 


N,(f) 


HD76932 


5770±150 ± 100 


-I.0it0.16 ±0.2 


2.06±0.16±0.2 


< 0.03 


Sg 


(d),(9) 


HD 110283 


56G0±150 ± 100 


-2.G±0.1G ±0.2 


2.13±0.1G ± 0.2 


< 0.03 


Sg 





Figure 1: Graph of ^Li and ^Li abundances in halo stars; triangles denote upper limits, 
filled circles correspond to actual detections of ^Li, and open circles to ^Li abundances. 
Error bars are 95% c.l., obtained by summing quadratically statistical and systematics 
error bars. Solid curves represent the standard evolutionary curve for ^Li and ''Li (no 
stellar component is included for ^Li), and the dashed curves delimit the 95% c.l. on the 
^Li prediction. 

Figure 2: Graph of ^Be abundances vs metallicity. The solid and dashed lines correspond 
to the ®Li evolutionary curve and its 95% c.l. limits (without any primordial component) 
shifted to the level of ^Be abundances. 

Figure 3: Histogram of the statistics of the ^Li survival fractions {g^) in HD84937 (heavy 
solid), HD160617 (light sohd), and HD201891 (dashed). 

Figure 4: Plot of ^Li abundances vs T^ff, after having corrected the observed abun- 
dances for their trend in metallicity, a slope 1 in our case, to the metallicity of HD84937, 
[Fe/H] =—2.1. The error bars on the abundances implied by this correction were not 
applied; other error bars are 95% c.l.. 
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Figure 1: Graph of ^Li and ^Li abundances in halo stars; triangles denote upper limits, 
filled circles correspond to actual detections of ^Li, and open circles to ^Li abundances. 
Error bars are 95% c.l., obtained by summing quadratically statistical and systematics 
error bars. Solid curves represent the standard evolutionary curve for ^Li and ^Li (no 
stellar component is included for ^Li), and the dashed curves delimit the 95% c.l. on the 
^Li prediction. 
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Figure 2: Graph of ^Be abundances vs metallicity. The solid and dashed lines correspond 
to the ^Li evolutionary curve and its 95% c.l. limits (without any primordial component) 
shifted to the level of ^Be abundances. 
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Figure 3: Histogram of the statistics of the ^Li survival fractions (ge) in HD84937 (heavy 
solid), HD160617 (light solid), and HD201891 (dashed). 
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Figure 4: Plot of ^Li abundances vs T^ff, after having corrected the observed abun- 
dances for their trend in metallicity, a slope 1 in our case, to the metallicity of HD84937, 
[Fe/H]=— 2.1. The error bars on the abundances implied by this correction were not 
applied; other error bars are 95% c.l.. 
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